On early onset of quark number density at zero temperature 
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We study a longstanding problem in lattice QCD at low temperature and nonzero quark chemical potential on 
an onset of the quark number density at // = m,/2. We introduce a physical parametrization of the eigenvalues 
in the reduction formula of the fermion determinant. It is shown that the parametrization reduces the quark 
number density operator to an expression with the Fermi distribution of the quark. For each configuration, the 
eigenvalues of the reduced matrix correspond to one-particle energy states of a quark. The gap of the eigenspec- 
trum of the reduced matrix corresponds to the gap of the energy states, which causes the //-independence of the 
fermion determinant for small // at T = 0. Once y, exceeds the gap, the quark number density becomes nonzero 
for each configuration, which causes the early onset of the quark number density. 
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Properties of QCD at finite temperature (T) and quark 
chemical potential (//) are central issues in particle, nuclear 
and astrophysics in order to understand various states of mat- 
ters in our universe. One important nature of QCD is the 
deconfinement transition. Recent simulations on finer lat- 
tices show that the deconfinement transition is crossover [j]], 
and occurs at T pc = 150 — 170 MeV depending on ob- 



servables |12|-|4|]. It is expected that the crossover extends to 
nonzero fi and turns into a first order phase transition at a crit- 
ical endpoint(CEP), which has been studied in many refer- 
ences, e.g. iSlfJ. 

The QCD phase diagram is expected to contain rich struc- 
tures also at low temperatures O, [Zl |8|] . In addition, the dis- 
covery of a pulser with twice solar mass calls the reliable 
equation of state of QCD. Properties of QCD at low T and 
nonzero // was investigated in lattice simulations with a Glas- 
gow method |I3-13], where they obtained a signal of a finite 
density phase transition. They also obtained an onset of the 
quark number density at // = m^/2, which disagrees with a 
phenomenological expectation that the baryon number density 
remains zero for //^ = 3/z < Mn. Although several attempts 
were made, e.g. lfl4l[l5ll . the solution is not yet obtained. In 
order to reveal the QCD phase diagram based on first principle 
lattice simulations, one must confront this problem. 

It was shown theoretically 1711 that the fermion deter- 
minant dot A(/x) is independent of /i at T = 0. Cohen showed 
the //-independence for isospin chemical potential case by 
considering the 70 times the Dirac operator. Adams showed 
it for quark chemical potential case by using a reduction for- 
mula of the fermion determinant. For lattice QCD, we showed 
the //-independence of det A (/it) by using the reduction for- 
mula iflill in two approaches. One approach was to evalu- 
ate the //- and T-dependence of a factor det A(/x) / det A(0). 
We found that the factor becomes insensitive to /x as T de- 
creases. The other one was to take the low temperature limit 
of det A(/x) with the aid of properties of the reduced matrix, 
a N t scaling law, the existence of a gap, pair nature and a re- 
lation between the gap and pion mass. We found that those 
properties makes det A(//) //-independent for /t < m^/2. 

Thus, the //-independence of the fermion determinant at 



low T and small // was explained both in theoretical stud- 
ies and lattice simulations. It is likely that the approaches 
in Refs. GUH| are closely related. Now, it is meaningful 
to consider underlying physics behind the early onset of the 
quark number density by extending the previous approaches. 

In this paper, we address the problem of the early onset 
of the quark number density by using the reduction formula. 
According to the previous studies, we introduce a physical 
parametrization of the eigenvalues of the reduced matrix Q. 
It will be shown that the parametrization reduces the num- 
ber density operator of the quark to an expression with the 
Fermi distribution. This makes it possible to study the prop- 
erty of the quark number density at T = in a familiar tech- 
nique. The gap of the eigenvalues of Q plays a similar role to 
a Fermi energy in the Fermi distribution. The /t-independence 
of det A(//) at T — for /x < m^/2 is explained by a simple 
argument based on the Fermi distribution with the energy gap. 
We will discuss the origin of the early onset. 

We employ clover-improved Wilson fermions of Nf = 2. 
The result can be applied to other actions if the reduction for- 
mula is available. First we divide the Wilson fermion matrix 
A into the spatial part and temporal part, 



A(/z) =B-k |e + "°(r - lA )U A {x)5 x , x+l 
+e-^{r + li )U\{x')5 xl ^ i , 



(1) 



i=i 

+ (r + li)U}{x')5 x , x _i - kCsw?>x,x> ^ Q'^uF^u, (2) 

where B is the spatial Wilson matrix, n the hopping parame- 
ter, r the Wilson parameter, Csw the clover coefficient. We 
introduce two block-matrices 



(3a) 
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>(x, y, U) rl v - 2k r^6(x - y) U?(y, U). 



(3b) 
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r± = (r ± 74,)/2 are projection operators in case of r = 1, 
which is used in the derivation of the formula. Using the block 
matrices, det A(/z) is given by IU9I1 



dot AO*) = C Nrcd/2 Co dot (£ + Q) , (4a) 

Q = Ufa 1 Pi), (4b) 
1=1 

N t 

Co = Y[det( ai ), 



(4c) 



where £ = e^/ T , TV = 4N c NfN t and AT red = 7V/iV t . N 
and iVped are the dimensions of A and Q, respectively. Q and 
Co are independent of /i. Calculating \Q — XI\ = 0, we obtain 



detA( A1 ) = C r Ar - d/2 n( A » + ^' 



(5) 



which is the analytic function of fj,, and provides values of 
det A(/*) for arbitrary \i. 

The reduced matrix Q has three important properties. One 
follows from the 75 hermiticity that the eigenvalues of Q ap- 
pear in pairs A and 1/A*, which can be seen in a two-peak 
behavior in Fig.Q] Accordingly, the eigenvalues are classified 
into two groups, |A| < 1 and |A| > 1. Here the numerical 
simulations were done in Nf — 2 clover-improved Wilson 
fermions and RG-improved gauge action with 8 3 x Nt lattice 



and with m ps /mv = 0.8. Data are taken from Ref. 01811 . 
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FIG. 1: The histogram of the eigenvalue distribution as a function 
of In |A|/iV t for Nt = 4(blue) and N = 8 (red) with a fixed scale 
a. The angular part is integrated out. The agreement between them 
implies the scaling law A ~ l Nt . The two peaks are a consequence of 
the pair nature A and 1/A*. No eigenvalue exists for |A| ~ 1 (gap). 

Second property is the gap in the eigenspectrum of Q. No 
eigenvalue exists in the vicinity of |A| = 1, see Fig. [2] Gibbs 
pointed out [20] that an eigenvalue near the gap is related to 
the pion mass. Later Fodor, Szabo and Toth investigated fur- 
ther the connection between the hadron spectrum and eigen- 
values of Q 112 111 . In Ref. 11811 . we have also calculated the 
pion mass in Nf = 2 Wilson fermion case. 

Third one is the Nt scaling law. Figure Q] shows the his- 
togram of A as a function of In | A|/iV f for Nt — 4 and N t — 8 
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FIG. 2: T-dependence of the largest eigenvalue in |A| < 1, i.e., 
rmiX|A„|<i |A n |, which is related to the pion mass and e m in. The 
ensemble average was taken over 400 configurations for mps /my = 
0.8 and 50 configurations for mps/mv = 0.6. 



with a fixed scale a. The agreement of the results indicates 
that an eigenvalue is parametrized as |A| ~ l Nt , where I is 
a real. In Ref. Ill 811 . we derived the low-T limit of det A(/i) 
by using these properties of A, and found that det A(/i) is fi- 
independent for fi < m^/2 in the low T limit. 

According to the N t scaling law, we introduce a 
parametrization of the eigenvalues of Q as |A| = cxp(— e/T) 
and 1 1/A* I = exp(e/T), where we describe larger half of the 
eigenvalues in terms of the smaller half of the eigenvalues. 
This parametrization was also implied in Ref. IU7I1 and sug- 
gested from an analogy between Q and the Polyakov line 11811 . 
The Polyakov loop describes the free energy of static quarks in 
heavy quark limit : (P) - er F / T with P = tiJJ^ (7 4 (^). 
Considering an eigenvalue VA, (| A < 1), the analogy between 
P and Q may suggest that A is related to an energy of a dy- 
namical quark. The eigenvalue A has the counterpart 1/A*. 
Because the pair nature results from the 75 hermiticity, it is 
natural to identify 1/A* as an energy of an anti-quark. 

Adams showed that the reduction formula reproduces the 
fermion determinant obtained from a Matsubara frequency 
summation method in free field case lfl7ll . We show this cor- 
respondence in the present case. Using the pair nature of the 
eigenvalues, Eq. (|5]l can be rewritten as 



iV red /2 



detA( A1 ) = C J] (A^l + Anr'Xl + A^)- (6) 

n=l 

For a while, we drop the phase of A. If we substitute A = 
e _£ / T and £ = e~^/ T , then Eq. (O is reduced to 

N Icd /2 

detA( M ) = C [] e e / T (l + e-( £ -^ T )(l + e -( e +^/ T ), 

n=l 

(7) 

which agrees with the Matsubara frequency summation 
method. Note that there is a degeneracy of Dirac compo- 
nents in the free field case. This correspondence suggests that 
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ea = — ln|A|/iV t describes one -particle energy states of a 
quark for a single configuration. 

Here, we comment on the ensemble average. The reduc- 
tion formula provides the fugacity expansion of the fermion 
determinant: det A(/i) = c n£, n , where c„ are functions of 
A„. Taking ensemble average correctly, a canonical partition 
function with a fixed quark number can be obtained, which 
provides a physical free energy of QCD. On the other hand, 
e corresponds to one-particle energy states of a quark for sin- 
gle configuration. This is perhaps related to a "quasi-energy " 
introduced in Ref. ifltill . which is eigenvalues of the 70 times 
Dirac operator. The ensemble average is crucial to calculate 
thermodynamical quantities, while it is convinient to consider 
e for the present purpose. 

Let us consider the number density operator. Equation (|6) 
is rewritten as 

/AW2 

detA(/i) =C exp [~ lnA » 

\ n=l 

+ ln(l + A„r 1 ) + ln(l + A;e)l I (8) 



For large volume, it is reduced to a spectral representation 

det A(/x) = C exp ( AN C V S f dXp(X) \- In A* 

V J w<i L 

+ ln(l + Ar 1 )+m(l + A*£)l ). (9) 



Here V s = N^, and p(X) is the spectral density defined on the 
complex A plane, p is normalized by 



Then, 



j d\p(X) = 1, 



' dXp(X) = / dXp(X) = 

|A|<1 J|A|>1 



(10a) 
(10b) 



Here the normalization conditions are satisfied in any N t . 
The number density operator for quarks is given by h 

[V^Tid/dfj.) det A(//)]/ det A{p). We obtain 



Ar 1 A*£ 



/ dXp(X) . 

|A|<i V 1 + A C 1 1 + A*C 



Now, we employ A = e £ / T + ie and £ = e ^/ T , 

1 1 



(11) 



11 = l dXp(X) 
|A|<1 



1 + e {e-ti)/T-%e I _|_ e (e+ M )/T+i9 J ' 

(12) 



and 



f dX= [°° %e~^ T P dS. (13a) 

J\\\<1 JO 1 J -IT 



de 

" = /o T e 



-2e/T 



1 



d8p(e-^ T+l9 ) 



1 



1 + e (e-n)/T-ie I + e (e+M)/T+ifi 



(14) 



where aside from 8, the first term in the large bracket is the 
Fermi distribution for a quark and the second term is that for 
an anti-quark. 

Once having the Fermi distribution, the zero temperature 
limit is obtained from a familiar technique. Taking Nt — > 00 
with a fixed lattice spacing a, the Fermi distribution is reduced 
to the step function. It is obvious from Nt -independence 
of Eq. ( fTOb that the factor e _2e / T /T does not cause any ill- 
behavior in Nt —> 00. Then, we obtain 



de 



-2e/T 



d6p(e 



e/T+i6\ 



(15) 



where p corresponds to a Fermi energy, and energy states be- 
low p are occupied at T = 0. 

The gap in the eigenspectrum of Q generates the minimum 
value e min given by maX| A |<i |A| = exp(-e min /T). The av- 
erage of e m in is shown in Fig. |2] No eigenvalue exists for 



p = for e < e„ 



Hence, we obtain 



n(p) = for p < e r 



(16) 



(17) 



Thus, the quark number density is zero at T = for small 
p for any configurations. According to Gibbs [20], e m i n 
corresponds to half the pion mass e m i n a = for 
each configuration. An alternative expression was found in 
Ref. 12 in . where the expression was given for the ensemble 
average. Although there is a difference between the results in 
Ref. B20I1 and 12111 . it is expected 12111 that the results become 
the same for large temporal lattice size. Hence, we employ 
£min = n is zero for configuration by configuration 

for p < e m in, and therefore its average is also zero, (h) = 0. 

Now, we discuss the onset of the quark number density. A 
phenomenological expectation is that the quark number den- 
sity starts to differ from zero at p = Mjy/3. This expectation 
is based on the fact that quarks are confined inside hadrons. 
The present result, e.g. Eq. (fl4] i indicates that the introduc- 
tion of the quark chemical potential actually affects an excita- 
tion of quarks even in the confinement phase at least for single 
configuration. For p < e m i n , p is not enough to excite quarks 
for any configurations, and the quark number density is zero. 
Once p exceeds e m j n , the chemical potential causes excitation 
of quarks and leads to nonzero quark number density for each 
configuration. This is the origin of the early onset of the quark 
number density. 

However, it does not necessarily mean that the average of 
the quark number density is nonzero. Once n takes nonzero 
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values, the phase of n also becomes nonzero. It is known that 
the sign problem becomes severe if \i goes beyond m w /2 II 1 811 . 
Hence, it is unclear what happens for m^/2 < /x < Mjv/3. 
The complex phase of n causes the cancellation of the quark 
number density over configurations, which may lead to (n) = 
for /x < M N /3. 

Another possibility is that finite density effects change e m i n . 
The relation e m ; n = m^/2 is based on the fact that the pion 
is the lightest hadron. Although QCD inequalities do not hold 
for finite density region, the pion is the lightest hadron at least 
for nuclear matter density. Hence the relation would hold at 
least for nuclear matter density. It seems that this possibility is 
unlikely. If (h) = for m^/2 < /x < Mjv/3, then it would be 
a consequence of the cancellation of the quark number density 
over configurations. 

In summary, we have addressed the longstanding problem 
of the onset of the quark number density. Based on previous 
studies, we considered it by using the properties of the reduc- 
tion formula obtained from the lattice QCD simulations. We 
parametrized the eigenvalues of the reduction formula based 
on some properties of the eigenvalues. The parametrization 
reduces the quark number density operator to the form with 
the Fermi distribution of a quark, where the eigenvalues of the 
reduced matrix correspond to one-particle energy states of a 
quark at single configuration level. Then, the /x-independence 
of the fermion determinant is explained by the familiar tech- 
nique. The present result shows that the quark chemical po- 
tential actually affects an excitation of quarks. Small /x can not 
excite any quark at T = 0, which causes the /i-independence 
of the fermion determinant at T = 0. On the other hand, 
if /x exceeds a certain value related to the pion mass, it can 
excite quarks even in the confinement phase at least one- 
configuration level, which is the origin of the early onset of 
the quark number density. If the quark number density re- 
mains zero up to [i = Mm /3, the complex phase of the quark 
number density would cause the cancellation over configura- 
tions. This result provides a hint to go beyond [i — ra^/2 in 
lattice simulations. 

For further confirmation, the spectral properties of the re- 
duced matrix should be investigated in future lattice simula- 
tions for smaller quark mass, and larger lattice. Particularly, 
the behavior of the gap and -dependence are important. It 
is also important to study the onset of the baryon number den- 
sity by both theoretical and experimental studies. Our result 
suggests the validity of the phase quench QCD and the orb- 
ifold equivalence at low temperature for /x < ra^/2 J22, 231. 
Phase quench simulations are also interesting for small [i. 
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